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1  Introduction 


Signal  representation  is  an  essential  step  in  many  engineering  signal  processing  ap¬ 
plications,  such  as  pattern  recognition  and  state  estimation.  The  key  task  of  signal 
representation  is  to  find  a  basis  to  parsimoniously  represent  a  given  signal.  In  ad¬ 
dition,  the  basis  should  (i)  capture  essential  signal  features  such  as  discontinuities, 
(ii)  match  the  smoothness  of  the  signal,  (iii)  accomodate  rapid  fluctuations  in  the  sig¬ 
nal,  (iv)  cater  to  the  stochasticity  and  the  distribution  of  the  signal,  and  (v)  be  usable 
on-line. 

If  the  signal  belongs  to  a  separable  space  such  as  the  Hilbert  space,  we  can  compactly 
represent  the  signal  by  a  countable  basis,  i.e.,  there  are  at  the  most  countable  infinity 
of  elements  (here,  the  basis  functions)  in  the  basis.  For  example,  suppose  the  signal 
emanates  from  a  second  order  stochastic  process,  i.e.,  a  process  whose  autocorrelation 
function  remains  finite  and  integrable  over  a  specified  interval.  The  space  of  second 
order  processes  — where  B  is  the  Borel-algebra  constructed  on  the  real  line 

R  and  p  is  the  underlying  probability  measure — is  complete  and  separable;  therefore 
we  may  employ  Karhunen-Loeve  (KL)  representation  to  obtain  a  basis. 

The  KL  representation  of  a  second  order  stochastic  process  y{t)  is  given  by 


y(0  = 


(1) 


where  Z  is  the  set  of  integers,  and  oij{t)  are  the  independent  solutions  of 


K(t,T)Qj{T)dT  =  Xjnjit) 


(2) 
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(i.e.,  aj{t)  are  the  eigenfunctions  of  the  covariance  function  A^(<,r)).  Since  A(t,r)  is 
self-adjoint  [9],  the  KL  representation  consists  of  expressing  a  given  signal  as  a  linear 
combination  of  a  few  orthonormal  modes.  Since  the  basis  is  orthonormal  and  consists 
of  the  eigenfunctions  of  the  stochastic  process  underlying  a  measured  signal,  the  KL 

representation  is  optimal  in  the  mean  square  sense. 

But,  computing  eigenfunctions  of  a  generic  second  order  process  is  extremely  te¬ 
dious.  In  addition,  from  an  implementation  standpoint,  the  KL  representation  needs 
a  few  ensembles  of  a  stochastic  sequence  to  estimate  K{t,  r).  In  the  absence  of  ensem¬ 
bles,  i.e.,  when  only  one  realization  is  available,  we  somehow  need  to  generate  some 
artificial  ensembles  in  order  to  develop  a  KL-like  representation. 

Alternatively,  Fourier  and  wavelet  representations  are  optimal  under  certain  condi¬ 
tions  and  are  commonly  used  in  practice  (3).  But  they  do  not  use  eigenfunction  basis. 
Of  all  the  non-eigenfunction  bases,  wavelet  bases  are  currently  the  most  attractive, 
because  of  their  ability  to  satisfy  the  five  desirable  properties  mentioned  earlier.  By 
proper  selection  of  wavelet  basis,  one  may  develop  optimal  representations  of  many 
classes  of  signals  [8,  11]. 

Wavelet  representation  is  algorithmically  simple  for  representing  transient  signals 
in  Aj-space.  Many  researchers  have  focussed  on  finding  optimal  wavelet  basis.  Coifman 
and  Wickerhauser  [12]  in  their  best  basis  formulation  proposed  that,  given  an  overcom¬ 
plete  set  of  basis  functions  belonging  to  a  wavelet  packet  (WP)  library,  we  may  choose 
their  optimal  combination  through  a  branch-and-bound  method  of  entropy  minimiza¬ 
tion.  However,  the  WP  library  and  hence  the  basis  functions  themselves  are  not  chosen 
optimally  with  respect  to  the  considered  class  of  signals. 
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For  the  pae.  one  decade,  aeveral  researchers  have  attempted  to  determine  the  op 
...  .avetet  has.  tor  d.erent  s.nat  Casses.  Ode.ard  et  at.  devetoped  opt.mat 
hnue  snpport  wavetet  basis  to  represent  hand-Umited  si.nats  hy  mtmmmms  t  e  ,n  n« 

„ormorthereso.„tio„erroroperatorateverysc.e»ithrespecttothescahn.^^^^^ 

of  that  scale.  Their  'robust”  basU  minimise  the  worst-ca^  error  tor  aU  s..n  . 
to  a  particular  trequeucy  band. 

unser  (111  showed  that,  when  the  signal  is  a  realisation  ot  a  stationary  process 

i.scjngtuncti„nistheimpu,seresponseotanide.ba„dpas^^ 

..ormiued  by  the  signal  energy.  His  analysis  showed  that  the  determination  ot 
optima  scaling  tunction  tor  hiorihoponel  ine.le.  represenlulmn  ot  a  generic  sign 

L2  space  is  mathematicaUy  intractable. 

Strintzis  1101  developed  the  necessary  conditions  for  a  globally  optimal 

..,o„otabm,d.bmitedmultidimensionalsignal,representedusingblorW  ' 

HO  tound  that  computing  a  glob.iy  optim.  scabng  tunctiou  tor  a  generic  multidimen- 
,onal  Signal  is  mathematic.,,  mtractabie  and,  in  many  cases,  not  re.icab,e. 

ComputationaHysimpiemethodsarecommonlyusedtodeterminesubopim- 

ingtunctionsUt,tO^Kither  the  structure  ot  the  synthesis  hi.  ^ 

pertorm  invert  wavelet  transtorm)  is  assumed  .  pnon,  and  the  analy 

,or  wavelet  transtorm,  is  constructed  theretromi  or  a  Butterworth.hlter.type  approa 

.oation  ot  the  ide.  ,inre.ir,ab,e  anaiysis-synthesis  pair  is  pertormcd. 

However,  most  ot  the  earlier  worhs  assume  the  Signal  to  emanate  troni  a  st^ 

.y  stochastic  process  with  a  hnown  variance.  None  ot  the  previous  worhs  ccp  ic,  y 
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inated  chaotic  signals  are  more  common  than  stationary  signals.  Hence,  research 
finding  the  optimal  representation  of  contaminated  chaotic  signals  seems  to  be  in 

der. 

In  this  paper,  we  develop  a  suboplimal  scaUng  function  using  the  local  eigentunc- 
tions,  extracted  as  described  in  the  toUowing  two  sections,  to  represent  signals  emanat¬ 
ing  from  a  chaotic  process  contaminated  with  low  leveU  of  noise.  In  other  words,  we 
develop  a  suboptimal  wavelet  packet  (WP)  representation  of  a  signal  which  is  a  single 
realization  of  a  chaotic  process  contaminated  with  low  intensity  noise. 

Our  methodology  consists  of  (i)  generating  artifical  ensembles  from  a  pseudo  probability 
space,  constructed  from  a  measured  signal,  as  explained  in  Section  2;  (ii)  extracting 
local  eigenfunctions  from  the  pseudo-probabability  space,  as  described  in  Section  3 
(these  local  eigenfunctions  may  be  used  to  develop  a  computationaUy  expensive  KL- 
like  representation  scheme-called  the  local  eigenfunction  representation-from  the 
obtained  ensembles);  and  (iii)  using  the  extracted  local  eigenfunctions  for  suboptimal 
WP  representation,  as  described  in  Section  4. 


2  Constructing  pseudo-probability  space  of  a 
chaotic  attractor 

Since  the  considered  signal  corresponds  to  a  single  trajectory  of  a  contaminated  chaotic 
process,  the  KL  representation  cannot  be  directly  used.  However,  by  appropriately 
constructing  a  pseudo-pwbability  sjxice  and  an  associated  proabability  measure  from 
the  TSD,  we  can  extract  artificial  ensembles  to  develop  a  KL -like  representation.  The 
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constructing  of  such  a  space  and  the  associated  probability  measure  for  a  dynamical 
system  is  studied  in  ergodic  theory  [4],  which  broadly  deals  with  developing  invariant 
measures  for  a  dynamic  system.  Before  we  proceed  further  towards  constructing  the 
pseudo-probability  space,  let  us  review  certain  relevant  aspects  of  dynamic  systems 
and  reinforce  the  concept  of  invariant  measures. 


2.1  Overview  of  ergodic  theory  and  our  approach 

Suppose  a  dynamic  system  is  modeled  in  terms  of  autonomous  [7]  nonlinear  stochastic 
differential/difference  equations  with  a  state  vector  x{t)  G  M  (the  state  space)  and  a 
discrete  output  y(n)  =  y{tn)  as  follows: 


dx  =  F(i)  dt  -b  £(x)  d£  ,  (^) 

y(n)  =  /i(i(n)) -|- v(n)  ;  n  =  0, 1, 2, . .  .TV  .  (4) 

Here,  (3)  represents  the  process  sytem  and  (4)  represents  the  measurement  setup.  Fur¬ 
thermore,  £(•)  is  a  nonlinear  stochastic  vector  field,  £(•)  and  h(-)  are  continuous 
transformations,  ^  is  a  vector  Wiener  process  [9]  that  accounts  for  dynamic  noise, 
i(n)  =  x{tn)y  and  v{n)  is  a  Gaussian  white  noise  sequence  (5]  that  accounts  for  ad¬ 
ditive  measurement  noise.  As  this  is  an  autonomous  system,  £(x)  and  £(x)  are  not 
explicit  functions  of  time  t.  Equation  (3)  is  an  Ito  equation  (5)  because  the  gam  matrix 
g{x)  is  a  stochastic  process. 

When  there  is  no  dynamic  noise,  i.e.,  £(x)  =  0,  the  solution  to  (3)  results  in  a 
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trajectory 


®(0  = 

where  /  ;  M  M  represents  the  flow  that  determines  the  evolution  of  *(<)  from  a 
specific  initial  condition  x(0).  If  the  sytem  is  dissipative,  i.e.,  the  volume  elements^ 
in  the  state  space  contract  as  the  system  evolves, ^  then  the  trajectory  generaUy  tends 
asymptotically  to  certain  compact  subsets.  If  S  C  M  is  a  compact  set  and  U  is  the 
largest  open  set  that  asymptotically  contracts  to  S  C  U,  then  S  is  called  an  attracting 
set  and  U  is  called  the  basin  of  attraction  of  S.  An  attracting  set  may  be  reduced  into 
certain  distinct  portions,  some  of  which  may  not  be  attracting. 

AU  disjoint  subsets  A,,  Aj,  A3, . .  •  of  S  that  are  attracting  are  caUed  attractors, 
while  S  -  (U>  Aj)  is  non-attracting.  An  attractor  A  C  M  is  associated  with  the 
following  four  properties  (7,  13]: 

1.  Invariance:  If  i(0)  G  A,  then  /(x(0),<)  G  A  Vf  >  0. 

2.  Attractivity:  For  some  U  C  M  defined  in  the  neighborhood  of  A,  £[U,  t]  C 

U  and  Urn, _oo /(«,<]  =  A. 

3.  Recurrence:  For  every  f  >  0  and  almost  every  i(0)  G  A,  3  <  >  0  such  that 

llx(O)  -  x(t)||  <  c,  in  effect,  the  trajectories  within  an  attractor  remain 
bounded. 

4.  Indecomposability:  An  attractor  is  disjoint  from  all  other  attractors  and 
'Volume  elements  refer  to  the  small  finite  chunks  of  state  space. 

-The  concept  of  dissipative  systems  emerges  from  an  Eulerian  perspective  rather  tlian  from  a  Lagrangian 
one. 
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tioned  properties. 

.Ue  flow  then  exhibits  a 

.t„ac.o,  A  may  ge.  separated  expeaeatially  as  the  system  evo  ^ 

known  as  the  sensiHee  dependence  on  initml  condition,  and  the  attrac  or 
nnued  a  s.mnye  nttmetor.  A  how  /C-)  -  o  parfmn.ar  initia,  condition,  ,s  sa,  o 

chaotic  it  the  trajectories  in  an  attractor  exhibit: 

sensitive  dependence  on  initial  conditions  bnt  are  bounded, 

2.  irregular  and  aperiodic  behavior,  and 

3.  continuous  broad  band  spectrum. 

of  dynamic  as  well  as  measurement  noise,  we  can  define 
Under  the  absence  of  dynamic 

thP  sLttractor  of  a  chaotic  process  as  follows; 
invariant  probability  measure  on  the 

■.u  r^rt  to  f(-  ■)  is  said  to  be  invanant  iff 
tku  eonlexl-de/ined  on  A  mU,  respect  ton.) 

„,/,A,.l)hoMs/or«n„tlcos.eoun.ch,ymon„un,ueso/.. 

It  the  trajectory  iies  in  the  attractor,  an  invariant  distribution  can  be  naturally  ob 

r.„od  because  the  attractor  is  a  compact  set  and  is  n, var, ant  under  the  no^  (^ 

r  attractor  is  evident  from  the  tollowmg 

The  existence  ot  an  invariant  measure  or 

theorem  [0]; 
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ApRned  on  a  compact  space  admits  at 
Theorem  1  Every  continuous  dynamic  system  defined 

least  one  invariant  probability  measure. 

as  The  probability  measure  may  be  extracted  usmg 

recurrence  theorem  stated  as  follows  [6]. 

,  Let  p-A^  A  be  a  measumble  transformation  on  a  probability  space 

Theorem  2  Let  p.  A  , 

.  a  Then  for  nlmost  all  £  €  ‘he  orkii 

(A,S„w).  Let  A  S  Ba  have  nM)  >  0. 

{p^{x))k>0  returns  to  A  infinitely  often. 


map 


£*'(x)  =  /U(0),tfc).«fc^0’ 


(6) 


.  U  ..e  U.C  a.  .....  ...e  cnU.  ..e  Po.nca..  .cUo„  A. 

by  Theorem  1,  ^  adn,i.s  a  measarc  .a  .^variant  «...  -pea.  to  P  . 

,  *  rr^nlps  of  the  probability  space  (A, 

demants  of  A  are  separate  iedependen.  sarr.p.e»  P 

where  6.  is  the  Bore,  algebra  constructed  on  A. 

„  ^  U  assumed  to  ..e  in  a  sma,,  neighborhood  of  £(0,  e  A  on  the  tratec 

I  e  A  =  cHB,(l(0))l,  a  Cosed  bai.  of  sine  c  constructed  about  £(0),  the 

the  eiernents  Of  A  are  independent  re— 


Thus,  ioca.  evolutions  from  neighborhood  points  of  a  single  trajectory  through  x(0) 
may  ' 


treated  as  indep . lent  random  trajectories  emerging  from  the  neighborhood  of 
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x(0).  In  addition,  for  a  chaotic  process,  the  independent  random  trajectories  emerging 
from  the  neighborhood  of  a  certain  initial  point  x(0)  €  A  remain  locally  close  together. 
Therefore,  (A,Sa,/^a)  ‘s  a  valid  probability  space  with  the  trajectory-segments  emerg¬ 
ing  from  A  as  realizations  of  a  stochastic  process  deftned  thereon. 

However,  when  the  dynamic  system  is  contaminated,  the  probability  space  con¬ 
structed  from  a  trajectory  thereof  is  not  valid  under  all  conditions  because  the  pres¬ 
ence  of  noise  may  destroy  the  structure  of  the  attractor  A.  Hence,  we  call  the  con¬ 
structed  probability  space  a  pseudo-probability  space.  Thus,  from  the  constructed 
pseudo-probability  space,  the  separate  realizations  consists  of  local  evolutions  from 
the  neighborhood  of  i(0),  all  belonging  to  a  single  trajectory. 

2.2  Pseudo-probability  space  construction  procedure 

In  the  foregoing,  we  assumed  knowledge  of  x(<)  €  A  C  M.  But  when  only  a  scalar 
TSD  y(u)  corresponding  to  the  process  output  is  available,  we  need  to  extract  samples 
using  the  concept  of  pseudo-probability  space,  which  in  turn  requires  the  topological 
notion  of  neighborhood.  For  every  y(n),  there  exists,  a  diffeomorphism  mapping  y(u) 
into  the  observable  portion  of  x{t).  Alternatively,  if  we  reconstruct  the  attractor  of  the 
original  dynamics  from  the  TSD  using  lag-coordinates  [2],  there  is  a  diffeomorphism 
connecting  the  reconstructed  state  vector  C(*)  observable  subspace  of  the  actual 
state  vector  x(*).  The  state  vector  ((•)  of  the  lag-reconstructed  dynamics  lies  in  a 
vector  space  of  dimension  df;,  the  embedding  dimension.  We  would  like  to  stress  that 
the  procedure  for  determining  an  estimate  of  dp:  and  obtaining  ^{n)  is  nontrivial.  For 
the  details  thereof,  the  reader  is  referred  to  (l). 
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Original  attractor 


Lag-reconstructed  attractor 


x2(n)  -20  -20  x1(n)  y(n+1)-20  -20  y(n) 


igure  1:  Rossler  attractor  obtained  from:  (a)  actual  solution  of  the  differential  equation, 
))  lag  coordinates. 

We  illustrate  the  construction  of  the  pseudo- proabability  space  through  the  follow¬ 
ing  example.  The  original  attractor  and  the  lag-reconstructed  attractor  of  the  Rossler 
attractor  corresponding  to  the  three  first  order  differential  equations 


iy 

-{X2  +  X2^) 

i  = 

X2 

= 

x\  +  0.15x2 

0.20  +  X3(xi  -  10.0) 

(7) 


with  y{n)  =  x\(tn)  ^  n  =  l,2,...iV  are  shown  in  Figure  6.  In  the  lag-reconstructed 
vector  space,  at  a  particular  point  ((n)  on  the  attractor,  the  set  of  points  ± 
£),  i  =  0,1,2,...}  resulting  from  the  evolution  (both  forward  as  well  as  backward) 
of  C(^)  constitutes  a  slmnd.  If  the  TSD  is  chaotic,  two  nearby  points  located  in 
different  strands  of  the  same  trajectory  (read  TSD),  locally  stay  close  together  before 
exponentially  diverging-off. 

Thus,  if  we  can  identify  a  set  of  neighbors  1  ^  ^  ^  ^h)  about  a 
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*  J.  of  the  reconstructed  trajectory, 
r^lfferent  strands  oi  vn 

,,ec>r.e.  po..  .  „„p..  ..e 

the  evolution,  of  these  neighbors  npp 

,,.,„„etions.  .he.engthonheen.mh.es..  e 

n  ehv  is  the  *cone.a.hm  length  Lot  the  TSDUl. 

Strand,  which  usually  procedure 

r  a  thus  we  may  obtain  the  eigenfunctions  using 

ensembles  generated  thus,  we 


as  in 


KL  representation. 


3  Extracting  local  eigenfunctions 

t  d  as  described  in  the  previous  section,  we  can 

From  the  nrtiftcal  ensembles  genera  e  ^„3embles  emerging 

.  u„is  functions  to  represent  L-leng  E 

„  ,,,,,,  loea.  eigentunc. 

ftom  A  .  BU(0)).  an  open-haU  ahou  .(  •  ^he  stepwise 

.fonstrom  the  artiheia.  ensembles  is  same  as  that  rnKLrp 

orethodology  of  eigenfunction  extraction  ,s  as  follows. 

,  .  ,m(x(0))|l  and  the  length  of  the  samples,  ..c 

1.  Fix  the  neighborhood  s,«  llB(x(0))l| 

the  strand  length  L. 

,  start  with  the  initial  point  g(0)  on  the  traiectory. 

,  oenerateartihclai  ensembles  as  described  in  section. 

..  .  nrtifical  ensembles.  The  procer 

4,  Extract  local  oigenfunrlions  from 

is  same  as  in  Kl.  ,..,.te»enlalion. 

•  ,  f>iul  carry  out  steps  (3)  and 

3.  NOW  shift  to  a  new  location  on  the  trarectory  and 

(1). 
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l«"e‘t>)-  wccaadeve\opal0C3l®'e'”'““'“°“ 

US, a*  *e  ^traced  local  e.sea  o  ^ 

■"  *  7';;;;,  „  we  ^0.  a  .o^deau,  .a.e  ae.bWbood 
,  .e  abclab  Ibe  balUo  a  sla..^^^ 

-eaaa>becor.caads.Oeecal,a.^^ 

oae  ele.ea.,,  H-e  -Ulas  reprea 

,.,.baa,la— 

,„„.«oas.  eatcacled  asla,  a 
p„,  cbaollc  slgaala,  »e  ,,,  la  aaba.aa.laled 

aafficieatly  large-daed  ball,  to  ref 

„Cbeccsaltso.aaa>erlcaUxpena,eatabt 


val  oerformance  evaluation 

3.1  Numeric  P  ,,0)  -  p.2012 , 2.6967 , 2.215cf. 

1W.  =  3  imtiallocationCW-l'* 

l.lbeaosalecatUaclol,wese.d  - 

we  loaad  Ne  -  «  aeareat  ae,g  oc  „„„ 

geaecated  10  easea,b,ea  tbece  ,oa 

ttegeaecatedeaac™b,ea.O„,,6e,geal»act,oa^ 

doadaaa.  aad  tbese  aloae  wece  used  loc  s.ga 

roariately  abiWag  Utese  6  doaaaaal  eig 

Byappropna  able  to  acca.ately  repteaeat  tb 

„g„al  as, ag  local  cabaeatWeaeatcb, 


I'icurc  :!■  Six  dominapit  local  eigenfunctions  for  the 
=  3,  lag  =  ^  and  L  =  168  (only  168  dalapoint-long 


llosslcr  attractor.  Here, 
eigenfunctions  shown). 


Nb 


40, 
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•  roc  that  the  local  eigenfunction 

.  V  «  of  tkc  total  signal  energy.  This  , mpl.es 

XSD«.thm8%of  The  accuracy  of  the  local  eigenfunct, on 

.  •  oflPouate  for  chaotic  signals.  The  accur  y 

representation  is  noise 

level  because  the  considered  signal  eman  is  the  computation 

„„„ever.onemalorara»hacl,onocalelgent«nc.,  P 

.i„e  reuulred  to  find  appropriate  sWfts  of  e,gent«nct,ons 


4  Subop.im.1  ™v....  PaeW  repre„....l<.» 

f  f‘.r^nQ  as  the  scaling  func- 
™W„ation  of  dominant  local  eigenfunct.o 

We  used  a  linear  comhma  from  a  chaotic  s.gnal  using  a 

«on.  suppose  rrruretfie  local  eigen, uncons  extract 

sufficiently  large  A  and 


1  -1. 
A 


(8) 


I'no-  fc*'*'  ensemble.  The  scaling 
m  •  nic;  of  n  (t)  while  representing 
where  are  the  coefficients  o  A  ) 

function  is  then  given  by 


jgD 


(9) 


resulting  function.  From  the  vector  unan  • 

,3.cm-c^h/.mc..onWmpse„do.prohahi,it,space. 
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Table  1:  Comparison  of  entropies  resulting  from  wavelet  packet  representation  of  a  signal 


Scaling  function 

Entropy 

Coiflet  1 

5.73 

Daubechies  4 

5.77 

Daubechies  20 

5.43 

Cosine  Packet 

r  2.63 

Suboptimal 

0.03 

We  used  the  extracted  scaling  functions  for  WP  representation  of  signals  from  the 
Rossler  attractor.  The  length  of  the  ensembles  and  hence  the  codebook  function  was 
set  to  L  =  168.  We  decimated  the  codebook  function  by  4  to  result  in  a  scaling 
function  42  datapoints  long.  The  computed  scaling  function  is  shown  in  Figure  6,  and 
a  comparison  of  entropt^  resulting  from  WP  representation  with  our  scaling  function 
against  the  entropies  resulting  from  WP  representation  with  other  standard  bases  [12] 
is  provided  in  Table  1.  The  entropy  values  of  suboptimal  representation  are  smaller 
than  those  of  other  standard  bases  by  an  order  of  magnitude.  This  implies  that  our 
suboptimal  basis  is  more  suited  other  standard  bases  for  representing  chaotic  signals. 

A  practical  validation  of  our  novel  representation  scheme  is  provided  in  the  following 
section. 


^Entropy  is  a  rneasure  of  parsimony  of  representation.  The  smaller  the  value  of  signal  entropy,  the 
IS  the  parsimony.  [12] 


greater 
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Normalized  response 
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Table  2:  Comparison  of  entropies  resulting  from  wavelet  packet  representation  of  a  signal 


Scaling  function 

Entropy 

CoiBet  1 

4.79 

Daubechies  4 

5.32 

Daubechies  20 

4.53 

Cosine  Packet 

3.84 

Suboptimal 

0.05 

5  Application  to  signals  from  machining  sen¬ 
sors 

We  used  the  suboptimal  wavelet  basis  to  represent  the  acoustic  emission  (AE)  signals 
in  machining.  For  a  partcular  signal,  with  L  =  76,  Nb  =  50  and  ds  =  9,  we  computed 
50  ensembles  as  shown  in  Figure  6,  and  determined  the  local  eigenfunctions  therefrom. 
Only  3  eigenfunctions,  shown  in  Figure  6,  were  dominant.  From  these  eigenfunctions, 
we  constructed  a  scaling  function  as  described  earlier  in  this  paper.  The  use  of  this 
scaling  function  substantially  reduced  the  entropy  of  representation  as  revealed  in  Ta¬ 
ble  2.  Clearly,  the  entropy  of  WP  representation  with  the  suboptimal  scaling  function  is 
about  a  magnitude  less  than  that  for  WP  representation  using  other  scaling  functions, 
implying  the  practical  effectiveness  of  our  method. 


6  Conclusion 


We  have  thus  developed  a  novel  representation  scheme  for  contaminated  chaotic  signals. 
It  is  probably  for  the  first  time  that  a  repre.sentation  scheme  for  contaminated  chaotic 
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AE  signal 


Figure  5:  Artificial  ensembles  generated  from  an  acoustic  emission  signal.  The  artificial 
ensembles  were  generated  with  Nb  =  50,  d/r  =  9,  lag  =  G  and  L  =  GO. 
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signals  has  been  developed.  This  parsimonious  suboptimal  representation  scheme  may 
be  applied  to  various  engineering  disciplines  where  chaotic  signals  occur,  and  specifi¬ 
cally  to  the  machining  process.  In  fact,  we  have  met  with  reasonable  success  in  using 
the  scaling  functions  developed  in  this  work  to  extract  signal  features  from  AE  signals 
for  tool  wear  estimation  [2]. 
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